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ABSTRACT 



Lattice QCD should allow a derivation of the AI = 1/2 rule from first principles, 
but numerical calculations to date have been plagued by a variety of problems. 
After a brief review of these problems, we present several new methods for calcu- 
lating K ^ nil amplitudes. These are designed for Wilson fermions, though they 
can be used also with staggered fermions. They all involve a non-perturbative 
determination of matching coefficients. We show how problems of operator mix- 
ing can be greatly reduced by using point-split hadronic currents, and how CP 
violating parts of the K ^ im amplitudes can be calculated by introducing a 
fake top quark. Many of the methods can also be applied to the calculation of 
two body non-leptonic i?-meson decays. 



1 Introduction 



One of the least well understood features of hadronic physics is the AJ = 1/2 rule in non- 
leptonic kaon decays. Decays in which isospin changes by A/ = 1/2 are greatly enhanced 
over those with A/ = 3/2. The particular example we focus on here is the ratio of amplitudes 
for K —>■ Tin decays: 

A{K ^ n7i[I = 2]) ^ ' ^ ' 

Although the origin of this large enhancement is not well understood, we do know that, in 
a QCD-based explanation, most of the enhancement must come from long distance, non- 
perturbative physics. This is because the contribution from scales where perturbation theory 
is reliable, say p ^2 GeV, is known to enhance the A J = 1/2 amplitude by only a factor of 
about two. Attempts to understand the remainder of the enhancement using models of non- 
perturbative physics have had partial success but we are far from having a convincing 
demonstration that QCD does explain the A J = 1/2 rule. 

In principle, lattice QCD is well suited to study this issue p[. Indeed, practical approaches 
have been developed for both Wilson 0, ^ and staggered fermions |^ . The problem is that 
these methods have not yet yielded useful results. Indirect methods, based on using K —>■ n 
and K —^vacuum amplitudes, have large statistical errors 0, H, while direct calculations 
of K ^ TTTT amplitudes have been done at unphysically large quark masses, for which the 
presence of a scalar resonance may distort the two pion signal |, |T^. There has been 
little work on the problem in recent years. 

In this paper we revisit the lattice approach for the case of Wilson fermions (or 0{a) im- 
proved versions thereof). We reevaluate the existing methods, and propose a variety of new 
approaches. These vary from a method of comparable simplicity to that of Bernard et al. Q , 
to more speculative ideas that are likely to require much smaller lattice spacings than those 
presently available. The latter are needed to calculate the CP violating parts of the K ^ titi 
amplitudes. We also reappraise the calculation of ^ vr matrix elements in the light of the 
recent progress made in non-perturbative renormalization techniques. 

Existing methods rely on lowest order chiral perturbation theory to connect the unphysical 
amplitudes which are calculated on the lattice to the desired physical amplitude. One ad- 
vantage of some of our new methods is that they are not dependent on chiral perturbation 
theory. This means that they can, in principle, also be used to study 5-meson decays. 

The outline of this paper is as follows. In the next section we summarise the problems faced 
by lattice calculations, and then, in sec. |, recall how these are avoided by the direct method 
of Bernard el al. With the stage thus set, in sec. ^ we then present the first three of our 
new methods, which also involve the calculation of K tttc amplitudes. This is followed 
in sec. |^ by a reappraisal of the indirect method using K ^ it amplitudes. In sec. ^, we 
present a more speculative method for studying the Al = 1/2 rule, which is based on the 
short distance expansion of the T-product of two weak currents. Finally in sec. |^ we suggest 
a possible idea for a non-perturbative determination of the CP- violating part of the AS* = 1 
weak Hamiltonian by introducing a fictitious top quark. Section ^ contains our conclusions. 
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2 The problem 



In this section we briefly review the source of the difficulty in calculating A{K vrvr) with 
Wilson fermions. Further details can be found in refs. |rT| and references therein. 

For scales below M^, but above the charm quark mass, the AS* = 1 part of the effective 
weak Hamiltonian can be written as 



njAS=l _ X 



0(±) 



C+(/i,M^)OW(/i) + C_(/x,M^)0(-)(/i)J , (2) 
(S7^f/)(M7» ± (s7»(M7^rf)] -[u^c], (3) 



where 7^ = 7^(1 — 75)72 and = VudV*^. Here and in the following we use the Euclidean 
metric. We are ignoring for the moment the contribution which arises when the top quark 
is integrated out. This is suppressed by Xt/Xu, where = VtdV^*, and, in the CP conserving 
sector, makes a very small contribution to the decay amplitudes. 

The operators O^"^^ have different transformation properties under isospin. In particular, 
0^~^ is pure I = 1/2, whereas O*-^-* contains parts having both J = 1/2 and I = 3/2. An 
explanation of the AJ = 1/2 rule thus requires that the K — > vrvr matrix element of C-0^~^ 
be substantially enhanced compared to that of C+O'-^-'. The short distance, perturbative 
part of the enhancement is that provided by the ratio of Wilson coefficients, C_/C+, while 
the long distance, non-perturbative contribution comes from the ratio of matrix elements of 
the operators. Consider first the short distance contribution. At ^ = M^, the two Wilson 
coefficients have nearly equal magnitudes, |C_/C+| = 1 + 0[as{M^)]. The enhancement 
arises from the renormalization group evolution down to /i ~ 2 GeV, at which scale one finds 
|C_/C+| ~ 2. This factor is too small by an order of magnitude to explain the A/ = 1/2 
rule. 

The remainder of the enhancement must come from the matrix elements of the operators, 
and these are the quantities that we wish to calculate on the lattice. In fig. we show the 
Wick contractions that contribute to such matrix elements. The only part of Ti-^^^ which 
gives rise to A J = 3/2 transitions is the 7 = 3/2 part of 0^~^\ and for this operator only 
contractions (A) and (C) contribute. In contrast, all four contractions are non-vanishing for 
the 1=1/2 parts of 0^'^\ Thus, in order to reproduce the AJ = 1/2 rule, the sum of the 
contributions of (B) and (D) must be an order of magnitude larger than those of (A) and 
(C). In fact, as discussed below, contractions (C) and (D) give non-leading contributions in 
chiral perturbation theory compared to (A) and (B). Thus we expect contraction (B) to be 
the "source" of the AJ = 1/2 rule. 

There are two major difficulties which arise in the calculation of these diagrams using lattice 
QCD. 

1. Decay amplitudes into two or more particles cannot be calculated directly in Euclidean 
space. This follows from the theorem of Maiani and Testa |T2[. One must, instead. 



use a model to analytically continue correlations functions from Euclidean to physical 
momenta. The only exception is for final state particles at rest relative to each other, 
in which case there is no phase generated by final state interactions. This problem 
afflicts both AJ = 1/2 and 3/2 amplitudes, although it is likely to be worse for the 



former since final state interactions are stronger for two pions having 1 = [O] 
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(A) (B) 




(C) (D) 

Figure 1: Contractions contributing to the AI = 1/2 amplitude. The lines represent prop- 
agators on a background gauge field, and the resulting contraction is to be averaged over 
gauge configurations with the appropriate measure. The dot represents the operators O^^^ 
Different colour contractions are not distinguished. 



2. The "penguin" diagrams (B) and (D) allow the J = 1/2 operators to mix with op- 
erators of lower dimension with coefficients which diverge as inverse powers of the 
lattice spacing. This mixing leads to contributions to the amplitudes which are lattice 
artifacts, and must be subtracted. One must also account for the mixing with other 
operators of dimension six, although this is less difficult because the mixing coefficients 
do not diverge in the continuum limit. 



In the remainder of this section we expand upon the latter problem. We consider only 
the negative parity parts of 0^'^\ since the parts with positive parity do not contribute to 
K ^ TTTT amplitudes. 

Lower dimension operators which can appear on the lattice must have the same flavour and 
CPS parity as O*-^-*. CPS is the transformation obtained by combining CP with (s ^ d) 
interchange |T^. Both O*^^-* have positive CPS parity. These symmetries allow mixing with 



only two lower dimension operators (aside from operators which vanish by the equations of 
motion): 

Op = (m, -md)s75rf, (4) 
Oa = go{ms-m,d)so^yG^ud, (5) 

Here g^ is the bare coupling constant and G^^ the dual field strength. The factors of nis—md 
are required by CPS symmetry. At dimension six, CPS forbids mixing with other four- 
fermion operators, and SU(4) flavour forbids mixing between 0^~^^ and 0^^\ It is noteworthy 
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that the allowed mixing for the negative parity parts of O^^^ is, up to this stage, exactly as 
in the continuum, despite the fact that chiral symmetry is explicitly broken on the lattice 
i- 

The mixing with Op and is further constrained by the GIM mechanism, which implies 
that the mixing coefficients vanish identically when ttIc = rrtu. In continuum perturbation 
theory, chiral symmetry requires that the mixing coefficients are quadratic functions of the 
quark masses and thus vanish as ml — rn^. On the lattice, by contrast, chiral symmetry is 
broken, and the GIM cancellation gives rise only to the factor ttIc — mu- This additional 
factor is sufficient, however, to make the mixing of 0(±) with an effect of 0(a). bmce m 
this work we are not attempting to remove 0{a) effects from the matrix elements of O'-^^ 
we do not need to consider mixing with Og-. We therefore arrive at the following form of the 
renormalized operator 



0(^)(/x) = zW(/za,^o') 



0^^\a) + {m,-mu)^Opia) 

ct 



+ 0(a). (6) 



HereO(±)(a) and Op{ a) are bare lattice operators, C^^-* are the mixing coefficients, and Z'-^-' 
are the renormalization constants which cancel the logarithmic divergence of the bare four- 
fermion operator. The precise definition of the lattice quark masses in eq. is unimportant, 
since in all practical methods the entire coefficient of Op is determined non-perturbatively. 

We now show that, in spite of the fact that it multiplies a linear divergence (see eq. @), 
it is sufficient to determine C^^-* with an error of 0(a), to obtain the physical amplitudes 
to the same precision. We first note that, in the continuum. Op is proportional to the 
divergence of the axial current, apart from terms which vanish by the equations of motion. 
Thus, in the continuum, it does not contribute to on-shell matrix elements for which the 
momentum inserted by the operator, Ap, vanishes. This is why the mixing of 0(±) with 
Op in the continuum is irrelevant to physical amplitudes. On the lattice, however, the 
breaking of chiral symmetry leads to 0(a) corrections to the PCAC equation even after 



renormalization ||15| 



{h\d^A^ - (m, + m,)P\h) = {h\XA\h2) = 0(a) . (7) 

Here = Z aS'') ^j^'^^d and P = Zpsj^d are the renormalized operators, ma and are 
renormalized quark masses, 9^ is a lattice derivative, and hi 2 represent hadronic states. 
Thus the K hit matrix element of the subtraction term is 

0(=^) _ O(^) ms-rrid I 

-^—{mc-mu){ms-md){7iiT\sj5d\K) = -^—{nic-mu)^—^ — {7!-]T\d^A^-XA\K) . (8) 

a a TUg + nid Zp 

The d^A^ term leads to the divergent contribution proportional to Ap/a. But when Ap = 
0, the term proportional to X^/a remains, and is of 0(1) up to logarithms. Thus the 
subtraction is still necessary, but the magnitude of the subtraction does not diverge in the 
continuum limit, and implies that it is sufficient to know C*^^-* with an error of 0(a). 

In summary, to calculate the physical K -kti matrix elements of O'-^^ one needs both a 
model to do the analytic continuation from Euclidean to Minkowski space, and a method to 
subtract the operator (m-c — mu)Op with appropriate coefficients O^^-'. Having summarised 
the problem, the remainder of the paper is devoted to illustrate several methods which can 
be used to extract the physical amplitudes. 
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3 Calculations with mg = irid. 



Bernard et al. proposed an ingenious solution to both problems [Q. They suggest working 
with rris = ma, and calculating the Euclidean amplitude in which all three particles are at 
rest: 



•^ms=md 



(7r(pi = 0)7r(p2 = 0)|OW(/i)|ir(pK = 0)) . (9) 

ms=ma 

Setting rUs = rrid causes Op to vanish identically, and so removes the need for the subtraction. 
Working with the two pions at rest solves the problem of final state interactions. The final 
step in the method is to extrapolate from the unphysical amplitudes J^^^ to the physical 
amplitudes ^phys using lowest order chiral perturbation theory. The result is that 



2 2 
"^X.phys ~ ^7r,phys a{±) 

n 2 •^ms=md 



1 + 



m 



K 



X 



(10) 



where mK,ia.t is the mass of the lattice "kaon" used in calculating A^^^^^^. Due to the choice 
= fTT'd, this is also the mass of the lattice "pion" . ^ 4:71 f^^ is the scale which determines 
the size of higher order terms in chiral perturbation theory. Such terms typically give 25% 
corrections if mK,\s,t ~ "^i^.phys- The direct method should thus be adequate to test the 
A/ = 1/2 rule at a semiquantitative level. 

The main advantage of this method is its simplicity. It determines the coefficient of the 
contribution to A^^^^ of leading order in chiral perturbation theory without the need for a 
subtraction. A further simplification is that one only needs to calculate the contractions 
(A) and (B) in fig. [|, because, due to CPS symmetry, the other two contractions vanish 
identically when = md- In fact, since one can show that the physical amplitude depends 
on quark masses only at non-leading order in chiral perturbation theory, the contributions 
of diagrams (C) and (D), which are proportional to — m^, must be of non- leading order. 

The main disadvantage of the method is its dependence on chiral perturbation theory. Such 
dependence is unavoidable once one has set = rrid- Of particular concern is the fact that 
the final state interaction phase in the / = two pion channel is substantial at the physical 
point s 



m 



K 



(where s is the square of the two-pion centre-of-mass energy), even though it 
is a non-leading effect in chiral perturbation theory. Thus it is possible that the corrections 
in eq. (|10D are larger than the estimate mj^/A'^ |jl3|. In principle one could reduce this error 
by using non-leading order chiral perturbation theory. This, however, requires knowledge of 
coefficients which are unavailable from experiment ||16i . 



4 Alternative methods 

We wish to develop methods which allow us to reduce, and ultimately remove, the dependence 
on chiral perturbation theory. To do this we must use quark masses which are closer to their 
physical values, and in particular we must consider rris > rrid- Working with non-degenerate 
quarks implies that the mixing with Op is present, and must either be subtracted, or shown to 
be unimportant. We have devised a number of new methods which accomplish this goal. In 
this section we describe, in order of increasing complexity, the three which require relatively 
small changes in methodology compared to the proposal of section |^. More speculative 
methods are described in subsequent sections. 
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Before describing our methods, we first comment on the the overall renormalization factors 
Z^^\lJLa, Qq). These must be calculated for all the methods described in this section and 
in sections ^ and ^. In the past Z^^'' have been determined using one-loop perturbation 
theory. Although this may be satisfactory for a test of the A J = 1/2 rule, a non-perturbative 
determination is clearly preferable. We wish to point out that such a determination can be 
made by a straightforward extension of the non-perturbative methods used to normalize 



AS* = 2 operators using quark states [^. The key observation is that the anomalous 
dimensions of O*^^-* are unchanged by mixing with the pseudoscalar density, and thus coincide 
with those of the operators 

= (?/'i7^?/'2 ^/'37m75^4 + ^i7m75V'2 ^37m^4) ±(2^4), (11) 

where the subscripts 1,2,3,4 label four different quark flavours. Only open current-current 
diagrams need be calculated; penguin-type diagrams with quark loops do not contribute. 
The non-perturbative calculation of the renormalized O*-^^ is currently underway |18 . 



4.1 Method 1 

We first describe the method and then explain why it is valid. The ingredients are 

• Work with a non-perturbatively 0{a) improved fermion action, for which there are no 
errors of 0{a) in the spectrum, and the on-shell amplitudes of the improved currents 
obey the continuum chiral Ward identities up to O(a^). This can be accomplished 
using the methods of ref. |]T^. 

• Use the lattice four-fermion operator, suitably normalized, but do not subtract the Op 
term. In other words, set C^^^ = 0. 

• Choose quark masses such that itlk = 2m.,^. 



Calculate the K —>■ amplitudes with all particles at rest. Given the choice of meson 

i(±) 



masses this means that Ap = 0. The resulting amplitudes are denoted A^"^^ 



Determine the physical amplitudes using lowest order chiral perturbation theory, which 
gives 



■^phys 2 2 '^m!f=2m^ 



2 

1 / I \ I 1 1 1 ^ ~ \ 

(12) 



1 + 



m 



K_ 



The validity of this method is based on the following observations. If one uses fully 0{a) 
improved fermions, and the fully 0(a) improved axial current, A^, and pseudoscalar density, 
, then the corrections to the PCAC equation (|^) are of O(a^) rather than of 0(a) |T^. 
Furthermore, the improved pseudoscalar density is simply proportional to the bare lattice 
operator = ZpS'-f^d. Thus the argument following eq. (0) now implies that the matrix 
element of the subtraction term is 

C{±) (7(±) _ I 

-^—{nic - mu){ms - md){-n'n\s-i^d\K) = -^—{rric - rriu)—^^ -—r{TX'n\d^Al\K) + 0(a) . 

a a rus + rud Zp ^ 

(13) 

The important point is that the discretization error on the r.h.s. of this result is now of 0(a) 
rather than of 0(1) [cf. eq. (H)]. This means that if we can set Ap = 0, so that the term 
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involving dfj^Aj^ vanishes, the subtraction leads to corrections only of 0(a). Since these are of 
the size that we are neglecting, it is not necessary to do the subtraction. Note that to make 
this argument, we do not need to know the form of A^, nor the value of the renormalization 
constants and Zp. 

Since we have chosen quark masses such that rriK = 2m7r, and chosen all particles to have 
p = 0, we do have Ap = 0. Furthermore, the use of pions at rest allows us to avoid final 
state interactions, as in the = rrid method. 

We view this new method as complementary to the method of section ^ It is only slightly 
more difficult to implement practically. One complication is the need to use fully 0(a) 
improved Wilson fermions. This is, however, becoming standard in numerical simulations, 
now that the necessary "clover" coefficient csw has been determined non-perturbatively 
for quenched QCD [|1^]. Note that the method of section ^ can be used equally well with 
improved fermions, but that this does not provide any particular advantage over unimproved 
fermions. The results for K —>■ tttt amplitudes will have errors of 0(a) in both cases. To 
remove these errors would require not only improving the action but also improving the 
operators O'-^^. 

A second potential technical complication is the need to include all four Wick contractions of 
fig. One can, however, consistently drop diagrams (C) and (D) since they are non-leading 
in chiral perturbation theory, while we are relying in eq. (|12D on the leading order term. 

A possible advantage of this method is that it requires a smaller extrapolation to the physical 
point, since one is using quark masses which are closer to those of the physical quarks. 
Furthermore, if one includes diagrams (C) and (D), one can perform a partial test of the 
convergence of the chiral expansion. Diagram (C) is straightforward to calculate, while the 
disconnected diagram (D) presents more technical difficulties. We expect, however, that 
since (D) is Zweig forbidden, its contribution will be smaller than that of (C). 



4.2 Method 2 

A potential problem with the previous method is that the momentum inserted will not, in 
practice, be exactly zero. If so, the part of the Op term proportional to c?^^^ will give a 
"small" but divergent contribution proportional to AE/a, where AE = mx — 2m^. We can 
mitigate this potential problem by determining O^^-* non-perturbatively. We propose to do 
so by applying the condition 

(0|OW(/x)|ir) = 0. (14) 

We then proceed as in the previous method, except now keeping the subtraction term. Note 
that implementation of this condition does not require knowledge of the overall normaliza- 
tions Z'(^). 

For this relatively small increase in effort we effectively determine the 0(1) part of 0^^\ 
and thus reduce the error from the incomplete subtraction of the Op term by one power of 
a. In other words, whereas in Method 1 ignoring the subtraction term leads to an error 
~ AE/a + 0(a), in Method 2 the residue after implementing eq. ([T^) is ~ AE + O(a^). 
Since there are other sources of 0(a) errors, the total error in the K ^ irn matrix element 
is ~ AE + 0(a). Thus the extrapolation to AE = will be considerably less delicate. 
Furthermore, since the method no longer relies on the improved PCAC relation, it can be 
used with unimproved, or tree-level improved, Wilson fermions. 
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We now explain how the condition (JT^) is justified. At leading order in chiral perturbation 
theory, the K ^vacuum matrix element takes the form (/^ ~ 132 MeV) 



(0|OW(/i)|K0) = . (15) 

The coefficients (J'-^-* do not contribute to the physical K ^ ttti matrix element (see eq. (P!7| ) 
below), and thus can be set to zero. A similar unphysical arbitrariness remains away from 
the chiral limit. The point is that we can always redefine the renormalized operators as 
follows |1^ 

0(±) — ,0^^^+F^^\mc-mu)Op. (16) 

where the finite coefficients F^^^ are of O(Aqcd) but otherwise arbitrary. This redefinition 
does not change the physical matrix elements of O*^^^ since Op is a total divergence. But it 
does affect matrix elements in which momentum is inserted, such as those appearing in the 
condition of eq. (|T^. This condition can thus be fulfilled by an appropriate choice of F*^^). 
In particular, it can also be used if the strange quark is replaced by the bottom quark. 

The possibility of redefining the renormalized operators by finite terms, as in eq. (|TBp, implies 
that dependence of the K ^ mi matrix elements on AE is of 0(1) and does not vanish with 
the lattice spacing as one would naively expect. 



4.3 Method 3 

The previous methods require the use of quark masses whose ratio ms/rrid differs from its 
physical value. To move closer to physical quark masses while holding AE = we must work 
with final states in which 7^ 0. This forces us to deal with the Maiani- Testa theorem. 
We propose doing so using the method of Ciuchini et al. (CFMS) [^]. Given an analytic 
parameterization of the two pion scattering amplitude (e.g. resonance dominance) CFMS 
show how, in principle, one can extract the magnitude and phase of the physical amplitude by 
studying the Euclidean amplitude as a function of time for a variety of final pion momenta. 

An assumption made by CFMS is the "smoothness" of the off-shell amplitudes. In the 
absence of resonances coupled to the final-state mesons, this requires that the K ^ irn 
amplitudes do not vary rapidly with Ap. If there is a nearby resonance, CFMS assume that 
it dominates the momentum dependence of the amplitudes. The "smoothness" hypothesis 
is, in this case, the assumption that the couplings of the resonance to the final-state parti- 
cles are smooth functions of the external momenta. It is then is possible to find a simple 
parametrization of the amplitudes which describes their rapid variation (due to the presence 
of the resonance) with Ap. 

Thus we propose using the method of CFMS after fixing C^^^ with the conditions eq. (p!^. 
This removes the terms proportional to Ap/ a which might violate the smoothness hypothesis. 
Since this method makes no use of chiral perturbation theory, it can be applied also B 
decays. The only restriction is that the GIM mechanism must be operative, which limits 
one to processes which do not involve top quark loops. 
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5 K ^ TT Matrix Elements 

In this section we reconsider the method suggested in ref. which uses the K n matrix 
elements of the positive parity part of the weak Hamihonian. This method rehes on chiral 
perturbation theory, and in this respect we expect it to be of comparable accuracy to the 
direct method of Bernard et al. and the first two of new methods introduced in the previous 
section. Since only single-particle states are involved, the advantage of the "K ir method" 
is that it is technically easier to extract the relevant matrix elements. The disadvantage is 
that the operator mixing problem is much more complicated than for the negative parity part 
of 7i^, making an accurate evaluation of the matrix element of the renormalized operator 
difficult. 

We begin by recalling how the physical amplitude is obtained from the knowledge of the 
properties of the K ^ tt amplitude exploiting the Soft-Pion Theorems (SPTs). At leading 
order in chiral perturbation theory the physical amplitude takes the form (for Ap = 0) 



fn 



(17) 



The coefficients 7^^-' appear also in the expression for the K ^ 71 matrix element 



(7r+(p)|OW(/i)|ir+(g)) = -5(±)^ + 7(^^p-g. (18) 

J IT 

By studying this matrix element as a function oi p ■ q one can, in principle, determine 7^^-*, 
from which we obtain the K tttt matrix elements up to corrections of order m|-/A^. 

The expression for the positive-parity components of the renormalized O^^^ is ||^ [| 



2\ n(±) 



sub 



ZW(/ia, 



i=l 



rric - niu) 



OJa) 



(19) 



Here O^^'' are four-fermion operators of dimension 6 which have different chirality from 0^^\ 



They are listed in refs. |Tl]], I21[]-[B3]. The remaining operators are of lower dimension 



Oo = sd . 



(20) 
(21) 



Comparing the result in eq. (|^) with that for the negative parity parts of 0*^^\ eq. (^, we 
see that CPS symmetry provides much weaker constraints for the parity conserving parts. 
In particular, the loss of the factor of — rrid means that the mixing with the scalar density 
diverges like 1/a^. 

^Although, for convenience, we use the same symbol, Z^^\ the values of the overall renormalization 
constant for the positive parity part is in general different from that for the parity violating part, as a result 
of the explicit chiral symmetry breaking on the lattice. The operators O'--^^ in this section also refer to their 
positive parity components. 
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Reference ||^ suggested the following approach for determining the mixing coefficients: use 
perturbative values for the coefficients of the dimension 6 operators and Oo-, but determine 
Cj^-* non-perturbatively. The latter determination is to be made by adjusting C^^-* until 
the momentum independent part of the K ^ tt matrix element in eq. (|l^) vanishes, i.e. 
6^"^^ = 0. This is the positive parity analogue of the condition ([T^, and can be justified by 
similar arguments. 

One problem with this approach is that a perturbative determination of mixing coefficients 
is often not reliable. This is exemplified by the case of AS* = 2 operators, where the 
operator with non-perturbatively determined coefficients has, to a very good approximation, 
the expected chiral behaviour, while that with perturbative coefficients does not |21, 23| . 



Nevertheless, this approach should be pursued, because it could provide semiquantitative 
results for the A/ = 1/2 rule. To date, no useful results have been obtained. 

There are three approaches which have been proposed to determine the subtraction coeffi- 
cients of four-fermion operators in a non-perturbative way: Gauge Invariant Ward Identities 



(GIWIs) 0, Ward Identities on quark states [23| and non-perturbative renormalization be- 



tween quark states at large external momenta PT| 



We start with a discussion of the GIWIs, reformulated in the light of the recent progress 



in the exploitation of Ward identities in the context of improved actions fl^]. As we shall 
explain below, this is the most promising approach in the presence of mixing with lower 
dimensional operators (and is more general than the present case). Consider the following 
Ward identity in the chiral limit: 

Y: {d.Aiiy) O,S(0) $(a;i, x^, ■ ■ ■ , = <f (xi, x^, ■ ■ ■ , x^)) , (22) 

yen 

where 7^ is a region of space-time containing the origin bounded by two hyperplanes 7/4 = —ta 
and 1/4 = tf, and / labels the flavour component of the axial transformation. $ represents 
a multilocal gauge-invariant operator, with xi, X2, ■ ■ ■ , all lying outside TZ. dO^'^^/da-^ 
denotes the variation of the operators under inflnitesimal axial rotations of the flelds. As 
shown in ref. 0, in the chiral limit, there is a unique choice of the coefficients of the operators 
which belong to different chiral representations, i.e. C*^^-*, C^^-* and C^^\ such that the 

subtracted operators O^^^ satisfy this Ward identity. By varying the points Xi,X2, ■ ■ ■ ,x„. 



eq. (|22D corresponds to an overdetermined set of linear inhomogeneous equations which, in 
principle, allow for the determination of all the mixing coefficients. Using the property that 
the axial current is conserved in the chiral limit, it is convenient to rewrite eq. (^2]) as follows: 



r/n(±)/'r|N 

J2{{Ai{y, t,) - A{{y, -Q) Oi^i{0) $(xi, x^, ■ • • , x„)) = -^{^^ Hx„X2, ■ ■ ■ , x^) • 
y 

(23) 

The above equation shows that there are no contact terms arising upon integration over 
y WTl . The absence of contact terms implies that we do not have to include the mixing with 



operators that vanish by the equations of motion. Although these operators do not contribute 
to physical matrix elements, in general (see below) they must be taken into account in the 
determination of the mixing coefficients. Notice that the coefficients determined in the chiral 
limit are sufficient to predict unambiguously the physical K -ktt amplitudes. 

The other two non-perturbative methods mentioned above use quark and gluon correlators 
in a fixed gauge, and either impose normalization conditions at large Euclidean momenta 
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pT| , or enforce the Ward Identities on quark Green functions |]2^. The problem with these 
methods is that they require the inclusion of two additional classes of operator^: 



Gauge invariant operators which vanish by the equation of motion. These do not 
contribute to on-shell matrix elements, but do contribute to the off-shell correlators 
used in these methods. In the present case there are three such operators of low enough 
dimension 



(24) 



The coefficient of the first is proportional to (mc — m^ja^ while that of the latter two 
are proportional to {jric — ruu). These operators also appear in the continuum, but 
they can be removed in perturbation theory by studying how the form factors behave 
as one goes on-shell. This is not possible in a numerical simulation. 

Operators which are not gauge invariant. These are, however, constrained to be either 
BRST invariant, or to vanish by the equations of motion ||25|. Note that there is 

There are two such 

(25) 
(26) 



an exact BRST symmetry on the lattice after gauge fixing p6 



operators of low enough dimension, and with positive CPS parity, 

^^(^ + md)d - s{-1p + msi^d 



[rrir 



00 + md)d - s(-^ + m,)^rf 



By varying the external momenta, and using suitable projectors [^, it may be possible, in 
principle, to separate the contributions of these operators. The inclusion of five additional 
operators, when the normalization conditions are imposed on quark states in momentum 
space, seems, however, to make the possibility of a non-perturbative determination of the 
mixing coefficients quite remote in practice. 

These problems can be avoided, by working in configuration space, rather than in momentum 
space. This can be achieved by studying (unamputated) quark Green functions of the form 



(27) 



where all the points are separated. In this case none of the operators which vanish by the 
equation of motion can appear. To enforce the Ward identities, one uses eq. ( ^3]) with a non- 
gauge invariant $(xi, ■ ■ ■ , X4) = ipi{xi)ip2{x2)i!z{.x^)'^i{.Xi)- The expectation value has to be 
evaluated in a fixed gauge. Although, in this case, we have the same number of coefficients 
to determine, we prefer the GIWI method, since it requires the evaluation of gauge invariant 
correlation functions only. 



6 Construction of 7-^^ from first principles 

In this section we propose a method which, in principle, avoids the difficulties caused by 
mixing with lower dimension operators, and which automatically gives the effective weak 

^ For a more detailed discussion of why such operators appear see ref. Q . 
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Hamiltonian with the correct normahzation. In addition, it allows one to construct an 
improved weak Hamiltonian, i.e. one having errors of O(a^), given only the improved versions 
of the weak currents. The method does not use chiral perturbation theory, and thus applies 
equally well to the AS* = 1, AC = 1 and Ai? = 1 parts of the weak Hamiltonian. The 
method is speculative in the sense that it is likely to require more computational power than 
is presently available, although we expect it to be practical with the advent of Teraflops 
machines. 

The standard construction of the non-leptonic weak Hamiltonian begins with the expression 

= 9wJ d'xDj,{x- M^)T [J,l{x)JIl{Q)] , (28) 

where 

^ w w 

is the Ty-boson propagator and J pi is the (left-handed) hadronic weak current. One then 
performs an operator product expansion (OPE) on the product of the two currents in eq. (^8|) , 
which is justified by the observation that the dominant contribution to the integral comes 
from distances <C For physical amplitudes, one obtains in this way 

{h\nZ\h') = ^EC^(/x,M^)M^-'^'(/^|0»(/x)|/^') , (30) 

where di is the dimension of the operator 0*^*^(/i), and functions Ci{n, M^) result from the 
integration of the Wilson expansion coefficients, Cj(x; //) (defined in eq. (|33D below), with the 
VT-propagator. Schematically, suppressing Lorentz indices, one has 

C,(/i,M^)M^-* = J d^xD^(x;M^)c,(x;/i) . (31) 

The 0^^\fi) are quark and/or gluon operators renormalized at the subtraction point yU. The 
functions Ci(/i, M^) are evaluated in perturbation theory and their running with /x is dictated 
by the renormalization group equation which follows from the /i-independence of the l.h.s. 
of eq. (13). 

The sum in the expansion (|30|) is over operators of increasing dimension. We consider in the 
following only operators with dimensions di < 6, since the contribution from operators with 
cij > 6 is suppressed by powers of 

All the intricacies of operator mixing in the definition of the finite and renormalized op- 
erators, 0^^\n), come about because the integrals in (|28| ) and ( [3l| ) are extended down to 
the region of extremely small x. The complicated mixing for the 0^^\jj)^s in terms of bare 
operators arises from contact terms when the separation of the two currents goes to zero 
(i.e. when |x| is of the order of a). The problem is particularly bad because chiral symmetry 
is broken by the lattice regularization. This observation suggests that a simple way to avoid 
these complications is to implicitly define the renormalized operators by enforcing the OPE 
on the lattice for distances much larger than the lattice spacing a. We imagine proceeding 
in the following way: 

1. Take the T-product of two properly normalized weak currents, JpL{x)Jpj;^{0). If required 
these currents can be improved. 



12 



2. Measure the hadronic matrix element {h\T[JpL{x)Jp];^{0)]\h') in a Monte Carlo simula- 
tion, as a function of a; for | si ^ in the region 

a<|x|<AQ^^. (32) 

3. Extract the numbers {h\0^''\n)\h') by fitting in the region (^) the x-behaviour of 
{h\T[JpLix)JlL{Q)]\h') to the formula 



IT 



JpL{x)JlM \h') = J2c,{x;f,){h\0^^\f,)\h'), (33) 



where the Wilson coefficients q(x; fi) are determined by continuum perturbation theory 
using any standard renormalization scheme. The scale fi should be chosen so that If/i 
lies in the range defined by eq. (|32D. Since we only consider operators of dimension 
6 or lower, the T-product differs from the right-hand side of eq. (|33D by terms of 
0(|xPAqqj-,), which is an estimate of the size of the systematic errors in this procedure. 
Note that in eq. (p3|) an average over the points x and —x (schematically J{x)J{0) — >• 
l/2( J(x) J(0) + J(— x) J(0))) is implied in order to eliminate from the OPE terms which 
do not appear in physical amplitudes because of the integration over x in eq. (^8]) . 
These terms, however, would appear on the r.h.s. of eq. ( P3| ) were we to not perform 
this average. 

4. Insert the numbers {h\0^^\jj,)\h') determined in this way into the expression for the 
matrix elements of Ti.^, finally obtaining eq. 



For the implementation of this procedure, what is required is the existence of a window, 
eq. (|32|), in which the distance between the two currents is small enough so that perturbation 
theory can be used to determine the expected form of the OPE, but large enough that lattice 
artifacts are small. These artifacts will be suppressed by powers of a/x. Clearly the existence 
of such a window requires that we have a sufficiently small lattice spacing. At the same time 
the physical volume of the lattice must be sufficiently large to allow the formation of hadrons. 

A few remarks may be useful at this point: 



• The method determines directly the "physical" matrix elements of the operators ap- 
pearing in the OPE of the two currents, i.e. the matrix elements of the finite, renormal- 
ized operators 0^*^(/i), without any reference to the magnitude of the W^-mass. Thus 
we do not need to probe distances of 0{1/M^) with lattice calculations. 

• If the action and currents are improved, then the resulting matrix elements of 0*^*^(/i), 
and thus of H^, will also be improved. 

• The //-dependence of the matrix elements of the operators 0^*^(/i) is given trivially 
by that of the (perturbative) Wilson coefficients, Ci{x;fi). It compensates the related 
/i-dependence of the functions Ci{fi, M^) in such a way that the l.h.s of eq. (^) is 
independent of the choice of the subtraction point. A similar comment holds for the 
dependence on renormalization scheme. 

• Unlike the methods discussed in previous sections, this approach automatically yields 
hadronic amplitudes that are properly normalized (in the renormalization scheme in 
which the Wilson coefficients appearing in eq. (|33| ) are computed). 
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We now discuss the feasibility of the method in more detaiL The critical step is fitting to 



the form predicted by the OPE, eq. (|33|) . Typically more than one operator contributes to 
the sum, so one must be able to separate the contributions using their different dependence 
on X. The operators of interest are of dimension 6, and thus have Wilson coefficients which 
vary logarithmically with x. At leading order the form is 

C2o_ 
as{l/x)Y^o ^ as ii) , , . 

) ^ 4^^° log(3;/^) + • • • ' (34) 

where 'Jq^ is the one-loop anomalous dimension of the operator 0^^\ and (3o is coefficient 
of the one- loop term in the /5-function. By contrast, the coefficients of lower dimension 
operators diverge as powers of 1/x (up to logarithmic corrections). Thus if lower dimension 
operators are present they will dominate at short distances, and it will be very difficult to 
pick out the matrix elements of the dimension 6 operators. If, on the other hand, only 
dimension 6 operators appear then it may be possible to separately determine their matrix 
elements. How feasible this is depends on how large a range of x we can use, and on the 
magnitude of the differences between the anomalous dimensions. 

Fortunately, in the cases of interest, there are no operators of dimension lower than 6. 
Consider, for example, the AS* = 1 part of Ti.^. The operators which can appear in the 
OPE are O^^-* [defined in eq. (|^)], and in addition 

O' = {ml - ml) rsiO; - X)^l:d . (35) 

The GIM mechanism requires O' to vanish when mc = m„, while chiral symmetry requires 
both the quarks to be left-handed and that the GIM factor be quadratic in the quark masses. 
Although this operator looks new, its negative parity part is, by the equations of motion, 
proportional to (m^ — m^) Op, where Op is defined in eq. (^), and its positive parity part is 
proportional to {ml — ml) (m^ + m^) Os- So these are the same operators we encountered 
in sections and except for the overall factors. Since O' has dimension 6,, its coefficient 
function depends only logarithmically on x. 

To determine the matrix elements using eq. (^) we need the anomalous dimensions, which 
for the three operators are 

7^+^=4, 7j"^ = -8, 7o = 16. (36) 



In fact, the contribution of O' to the r.h.s. of eq. (|33| ) can be determined separately since its 
matrix element does not require any subtraction and can be calculated directly. As for O^, 
since their anomalous dimensions are well separated from one another, it may be possible 
to extract the corresponding matrix elements and then construct the physical amplitude of 

An important element of the procedure proposed in this section is that since it is the contin- 
uum OPE which determines the operators which appear, these are restricted by continuum 
symmetries. This is because, for \x\ ^ a, the lattice OPE matches that of the continuum 
with discretization errors suppressed by powers of a/x. 

The previous discussion shows how we can, in principle, remove the contribution of the 
scalar and pseudoscalar densities from any hadronic matrix element. Not only does the 
method work for both positive and negative parity parts of H^, but it also works if the 
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weak Hamiltonian carries momentum Ap. Thus the simplest way to test the method may be 
to calculate K ^ n matrix elements and then use chiral perturbation theory to relate these 
to the physical amplitude, as described in the previous section. The problems of operator 
mixing described in sec. ^ do not apply to the new method. 

Of course, for physical matrix elements one does not need to worry about the subtraction 
of O'. This is because the matrix elements determined by this method are those of the 
continuum operator, up to discretization errors. 

We end this section with an observation on the computational feasibility of this approach. 
The main difficulty is to have a sufficiently large range of values of |x| in order to separate 
the contributions from the different operators, and yet to satisfy the condition (|32|). These 
constraints make the method difficult at present, and it will only be fully exploited when 
Terafiops machines become available. 



7 A propagating top quark 

For CP violating processes in kaon decays, or for B decays where top-penguin diagrams 
enter at a Cabibbo-allowed level, the strategies described in sec. H and || for the negative 
parity operators fail, because the GIM mechanism is not operative. In particular O^^^ mix 
with all the penguin operators (see below). This makes the calculation of the mixing matrix 
of comparable difficulty to that for the positive parity operators described in sec. ^. For 
positive parity operators, the analysis of section ^ still applies. The difference is that the 
mixing coefficients of the magnetic operator and scalar density become more divergent 0] 
and this can make the numerical determination of the renormalized operators less precise. 
In order to circumvent these problems we propose two methods involving a fictitious top 
quark (with mass m^) which is light enough to propagate on the lattice. 

The basic idea is to work with two different scales: the first, /i, is larger than nit, so that 
the corresponding operator basis is as in the previous sections (0*^^^); the second, /i', is 
smaller than rrit so that a full set of penguin operators is generated. The matrix elements of 
the operators renormalized at the scale /i are computed numerically following the strategies 
explained in sees. By matching the result to the amplitude expressed in terms of 

operators renormalized at /i', we extract their matrix elements. In this way, at least in 
principle, we can obtain the matrix elements of the penguin operators without directly 
computing them. We now present the details of this procedure. 

At scales /i' below m^, when GIM is not operative, the form of the AS = 1 effective Hamil- 
tonian is 



H 



AS=1 ^ ^ 



- XtC{fi',M^,mt)-Q{fi')] (37) 

where = VqdV*^, q = u,c,t (Ac is eliminated by using the unitarity relation Ac = — A^ — A^). 
Here Q contains the "penguin" operators 

Q(/i') = (Qi(/x'),Q2(/i'),---,Q6(/^')) (38) 
and C are the corresponding coefficients 

C{fi') = (Ci(/x', MJ, C^ifi', MJ, Csifi', M^,mt)..., C^ifi', M^,mt)) . (39) 
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A convenient basis of operators when QCD corrections are taken into account is [p3| - 



P6| 



Qi = {sd)(^v-A){cc)(v-A) , 
Q2 = (sc)(v_A)(crf)(V-A) , 
Q3,5 = {sd)(^V-A)J2(^(l)iyTA) (40) 

Qi = J2(^l)(v-A){qd)(v-A) 

Qe = -2^(sg)(s+p)(gii)(s-P) . 
q 

Qi'* and Q2'* ^^^6 the analogous operators to Qi and Q2 with the up- and top-quark replacing 
the charmed one. Here the subscripts {V ± A) and {S ± P) indicate the chiral structures, 
and the sum over quarks q runs over the active flavours at the scale fi'. For simplicity we 
ignore electroweak penguin and magnetic operators; it is straightforward to generalize the 
following discussion to include them. The coefficient functions appearing above have been 
calculated up to non-leading order in perturbation theory in refs. [^]-|^. The effective 



Hamiltonian relevant for AB = 1 decays is simply obtained by replacing the s quark with 
the b quark. 

The part of H^^^ proportional to A„ is the same as that considered above in eq. and 
has been the focus of discussion for much of the paper. We have simply re-expressed it here 
in the new operator basis, in terms of which 

0± = {Q-,-Q,)±{Q^,-Q2), (41) 
C±(/i,M^) = ^[C,ifx,Mj±C2ifx,Mj]. (42) 

This part of n^i=^ gives the dominant contribution to CP conserving K ^ tttt amplitudes. 
These amplitudes can be calculated using the methods presented in sees. 

The difficulties arise for the part of Ti,^^^ proportional to Xt, which gives rise to CP violation 
in kaon decays. This part contains the penguin operators Qi, whose matrix elements are 
not protected by the GIM mechanism. To write a renormalized version of these operators 
requires subtracting Op and 0^ with appropriate coefficients, and accounting for mixing with 
all the other operators Qj. This is true for both positive and negative parity sectors. The 
methods of this section are designed to avoid these difficulties. 

We do so by introducing a dynamical top quark so as to keep the GIM mechanism operative. 
This will, however, be a fictitious top quark with mass satisfying 

l/a :$>mt:>mc:> Aqcd • (43) 

In other words, our top is light enough to propagate on the lattice, but, hke the physical 
top, it is heavier than the charm quark. As explained below, by using the fictitious top we 
can extract the matrix elements {h\Qi{fi')\h') which can then be inserted into the expression 
for H^^^, eq. (|37|). In this respect the method is similar to that of sec. For purposes of 
illustration we will restrict the discussion below to negative parity operators. 

We begin with the two matrix elements 

M(/i,m,)^(/i|g,(/i)-g*(/x)|/i'), t = l,2 (44) 
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evaluated at a renormalization scale satisfying ^ ^ ^ rht. Since GIM is operative, the 
analysis of sec. |^ applies (with rrit). Thus we can define Qi{n) — Q\{n) in terms of 

bare lattice operators by 



Qj{a) - g*(a) + (m, - mt)^0.,{a) 



a 



(45) 



where i,j = 1,2 and the subtraction coefficients C^^^ are determined by enforcing 



CO) 

{0\Qjia) - Q]{a) + (me - mt)^Op{a)\K) = . (46) 

The Zij are related by a simple change of basis to the Z^^) of eq. (P), and can be calculated 
either perturbatively or non-perturbatively. In this way we can obtain A^j(/i, m^) from a 
lattice calculation, as a function of m^, for some choice of /i. A simple choice is yU ~ 1/a. 

On the other hand, we can also consider the same matrix elements for a renormalization 
scale rfit ^ ^' ^ m^. In this case, the GIM mechanism is not operative, and the matrix 
elements can be expressed in terms of the six operators which appear in eq. ( ^0|) 



A^,(/x,m,) = Zr^{^^\^i,m^){h\Q,{^,')\h') . (47) 
i=i,6 

The rectangular matrix Z^^ can be calculated perturbatively by matching the theory with 
and without the fictitious top quark. This is exactly the method used to calculate the 
coefficients C in the theory with the physical top quark mass, except that in the physical 
case one must simultaneously integrate out both the W boson and the top quark. Here we 
are effectively integrating out the W first and then the fictitious top quark. The results for 
Z^^ at non- leading order can be reconstructed from those computed in refs. |^^-|3^. The 



running in the 1-2 submatrix is particularly simple (since the operators Qs-Qq do not feed 
back into Qi and Q2) and one can show that 

Zij{^'a,go)= J2 Zik{^\^)Zkj{^a,gl) , (48) 

fc=l,2 

with (i, j) = 1, 2. Here Zik is the inverse of the 2x2 sub-block of the mixing matrix appearing 
in eq. (^). This submatrix does not depend on nit, but only on /i' and /i. 

We now choose a value of fi' and vary nit. The six matrix elements of interest, {h\Qj{fi')\h'), 
are obtained by fitting the right hand side of eq. (^) to J^i{^,mt) computed numerically as 



in eq. (|4^), and using the renormalization matrix calculated perturbatively. Since the 
dependence on is logarithmic, this will not be easy. The procedure is analogous to our 
use of the x-dependence in sec. ^ to separate the renormalized matrix elements of operators 
appearing in the weak Hamiltonian. Having determined these renormalized matrix elements, 
we can insert them into the expression (^7p for the effective Hamiltonian. At this point the 
constraint M^^ ^ can be removed since the Wilson coefficients of the operators appearing 
in Ti^^p^ can be computed perturbatively for arbitrary values of , including mt = mt. 

Before discussing the errors involved in this procedure, we make the following observation. 
In the effective theory where the top quark has been removed, provided that /x' ^ Aqcd, 
we can evolve the renormalized operators from one scale to another using perturbation 
theory. In particular we can obtain {h\Q j{ii)\h') from {h\Qj{ii')\h') for all six operators 
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using the 6x6 anomalous dimension matrix computed in perturbation theory with the 
effective Hamihonian where the top quark has been integrated out. Thus we can directly 
extract the matrix elements of the operators from Aii{fi,mt). Note, however, that 

an accurate determination of the matrix elements of Q{fi) (or of Q{fi')) requires that the 
typical scale, Ahh', of masses and external momenta appearing in the physical process h' ^ h 
be much smaller than nit. This is because in the matching procedure we neglect terms of 
0{Ahh'/mt). 

An alternative method, in the same spirit as the approach followed in sec. |^, is the following. 
We can avoid the need for any non-perturbative subtraction by separating the two currents 
and using a fictitious propagating top quark. Thus we directly match (/i|T( Jpi(x) jj^(O)) |/i')top, 
where the subscript indicates the presence of the fictitious top, to the formula 

{h\T{J,L{x)JlLm\h%,, = E c,{x■,^^',m){h\0'^H^^')\h') , (49) 



1=1,6 



where 



Jpl{x)JIl{0) = s(x)7p(l - 75)t(x)t(0)7p(l - l5)d{0) 

- s(x)7p(l - 75)c(x)c(0)7p(l - 75)rf(0) • (50) 

The coefficients Ci^2{x; fi' ,mt) = Ci^2{x;^') are the same as those in sec. ||. The coefficients 
Ci{x] fi' ,rnt), with i = 3-6 are complicated functions of the anomalous dimension matrix 



which can be worked out from the results of refs. p9| and and computed numerically. 



Both methods proposed in this section require small enough lattice spacings to accommodate 
a number of scales. Like the method of sec. ^ their full implementation is likely to require 
the next generation of supercomputers. 



8 Conclusion 

In this paper we have suggested a number of new approaches with which to study the 
A/ = 1/2 rule using Wilson-like fermions. These methods can also be used for staggered 
fermions. 

In order to obtain the physical K —>■ nir amplitude without relying on chiral perturbation 
theory, or to study decays such as -B — vrvr, one must learn how to extract information on 



final state interactions from Euclidean amplitudes. The method of ref. |^ might make this 
possible, but detailed numerical studies are needed to assess whether it is practical. 

The calculation of the CP violating part of K rnr amplitudes with Wilson quarks is very 
difficult. A completely nonperturbative method may require the addition of a fictitious top 
quark. 

We have also reevaluated the method of ref. 0] involving K ^ n amplitudes. This approach 
is likely to be more difficult because of the large number of mixing coefficients which have 
to be determined non-perturbatively. It may however provide complementary information 
to the results obtained with the K ^ mr method, and a check of the accuracy of chiral 
relations. 

Only numerical studies will be able to confirm or refute our present intuition that the K — >■ 
Tin methods are likely to provide the better results in the near future. 
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